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(Dated: August 8, 2008) 

We report a search for the standard model Higgs boson in the missing energy and acoplanar 6-jet 
topology, using an integrated luminosity of 0.93 fb _1 recorded by the DO detector at the Fermilab 
Tevatron pp Collider. The analysis includes signal contributions from pp — > ZH — ► vvbb, as well as 
from WH production in which the charged lepton from the W boson decay is undetected. Neural 
networks are used to separate signal from background. In the absence of a signal, we set limits on 
cr(pp — > VH) x B(H — ► 66) at the 95% C.L. of 2.6 - 2.3 pb, for Higgs boson masses in the range 
105 - 135 GeV, where V = W,Z. The corresponding expected limits range from 2.8 pb - 2.0 pb. 
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PACS numbers: 13.85.Qk, 13.85.Ni, 13.85.Rm 

The Higgs mechanism, postulated to explain elec- 
troweak symmetry breaking, predicts the existence of the 
Higgs boson, which has yet to be found. The CERN e + e~ 
Collider experiments placed a lower limit on its mass of 
114.4 GeV at 95% C.L. Q. Global fits to precision elec- 
troweak data suggest a mass of Mjj < 160 GeV at 95% 
C.L. 0. In this range, the Fermilab Tevatron pp Col- 
lider has significant discovery potential. Searches in the 
missing energy and acoplanar 6-jet channel have been 
published by CDF [3] and DO 0. This channel is sen- 
sitive to ZH associated production when the Z decays 
to neutrinos and WH production when the charged lep- 
ton from the W decay is undetected. The result in this 
Letter supercedes our previous work. As well as bene- 
fitting from more data this analysis uses artificial neural 
networks (NN) for heavy flavor tagging (b tagging) and 
in event selection. 

The DO detector is described in Ref. [5J]. Dedicated 
triggers selected events with acoplanar jets and large im- 
balance in transverse momentum, ($t), as defined by 
energy deposited in the DO calorimeters. After impos- 
ing data quality requirements the data correspond to an 
integrated luminosity of 0.93 fb _1 [f|. Time-dependent 
adjustments have been made to the trigger requirements 
to compensate for the increasing peak instantaneous lu- 
minosity of the Tevatron. The selection criteria therefore 
varied somewhat, but typically required > 30 GeV 
(where Ifl? is the imbalance in transverse momentum 
calculated using only well reconstructed jets) for jets 
reconstructed at the highest level trigger, and an az- 
imuthal angle between the two leading (highest px) jets 
of A^jehJeh) < 170°. 

Event selection requires at least two jets with px > 
20 GeV, < 1.1 (central calorimeter) or 1.4 < \rj\ < 2.5 
(end calorimeters) and A(f)(jeti, jet 2 ) < 165°, where r\ is 
the pseudorapidity measured from the center of the de- 
tector (rj = — ln(tan|), where 6 is the angle relative to 
the beam axis). Reconstructed jets are corrected based 
on the expected calorimeter response, energy lost due to 
showering out of the jet cone, and energy deposited in the 
jet cone not associated with the jet [7[. We require the 
distance along the beam axis of the primary vertex from 
the center of the detector (zpy) to be less than 35 cm, 
and at least three tracks attached to the primary vertex 
to ensure b tagging capability. We also require $t> 50 
GeV and Ht < 240 GeV (where Ht is the scalar sum 
of the transverse momenta calculated using only well re- 
constructed jets) to reduce the contribution from it back- 
ground. A significant proportion of W / Z+jets events in 
which the bosons decay into charged leptons are rejected 
by vetoing isolated leptons (electrons or muons). 

Signal samples of ZH — > vvbb and WH — ► tvgbb 
(I = e,;U, t) were generated for 105 < Mb < 135 GeV 



using PYTHIA [8|. There are two types of backgrounds: 
physical processes modelled by Monte Carlo (MC) gener- 
ators and instrumental background predicted from data. 
ALPGEN Q was used to simulate it production with up 
to four jets. Samples of IF+jets (W decays to all three 
lepton pairs for light jets jj, bb and cc jets) and Z+jets 
(including Z — > vv and Z — > t + t~ processes for jj, bb 
and cc jets) were also generated separately using ALPGEN. 
Diboson processes (WW, WZ and ZZ) were generated 
with pythia. The samples generated with ALPGEN were 
processed through pythia for showering and hadroniza- 
tion. Ncxt-to-leading order (NLO) cross sections were 
used for normalizing all processes (NNLO for ti). All 
samples were processed through the DO detector simu- 
lation and the reconstruction software. The trigger re- 
quirements were modeled using a parametrized trigger 
simulation determined from data. 

As b tagging is applied later, jets are required to be 
"taggable" , i.e. satisfy certain minimal tracking and ver- 
texing criteria; a jet must have at least two tracks, one 
with pt > 1 and the other with > 0.5 GeV, each with > 2 
hits in the silicon vertex detector, and ATZ(track, jet) < 

0.5, where A1Z = \J (A<p) 2 + (A?y) 2 , with (f> being the 
azimuthal angle. The fraction of taggable jets was inves- 
tigated as a function of pr, rj and zpy using a VF+jets 
data sample. Jets in the simulation are corrected by the 
ratio of taggabilities measured in data and in MC, which 
is found to depend only on -q. Correction factors of 0.97 
± 0.01 and 0.95 ± 0.03 (statistical errors) are used for 
the central and end calorimeters, respectively. 

For events originating from hard processes with gen- 
uine missing transverse energy, the 1/Lt, $t and 
(where T/t is the negative of the vector sum of the pr 
of all tracks) point in the same direction and are cor- 
related. However, dijet events in which one of the jets 
has been mismeasured typically have $t pointing along 
the direction of one of the jets. Instrumental effects pro- 
duce events that tend to have T£t and T/t misaligned. To 
reduce instrumental background, we require: 

• minlA^Tjet,)} > 0.15, 

where min{A</>i(|?T JetJ} is the minimum of the 
difference in azimuthal angle between the direction 
of IpiT and any of the jets; 

• $ T (GeV) > -40 x min{A^(^ T , jetj} + 80; 

• A0(^ t ,?t) < tt/2, where A^ t ,7t) is the differ- 
ence in azimuthal angle between the directions of 
$t and Tf<r ; 

• -0.1 < A{$t,Wt) < 0.2, where A($t^t) = 
[J/jT — 1/[t)/{$t + Wt) is the asymmetry between 
TpT and fix ■ 
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FIG. 1: A($t,]/It) for data, MC physics background and 
instrumental background in the signal region, before im- 
plementing b tagging. The final selection corresponds to 
-0.K A(^ T ,ft T ) < 0.2. 

The residual contribution of the instrumental back- 
ground is determined from distributions in A($tiWt) 
and A4>($T:7t)- The instrumental background peaks at 
A{$t,Wt) < because it is dominated by poor quality 
jets that are taken into account when calculating $t but 
not ]^t- Signal and sideband regions are defined as hav- 
ing A0(^ T ,y T ) < tt/2 and A0(# T ,y T ) > tt/2, respec- 
tively. The shape of the backgrounds from simulated pro- 
cesses, for both regions, are taken directly from the MC. 
We fit a sixth-order polynomial to the A($tiWt) dis- 
tribution in the sideband region to determine the shape 
(before 6-jet tagging) for the instrumental background 
(after subtracting the MC background contribution) and 
a triple Gaussian for the signal region. We then do a com- 
bined physics and instrumental backgrounds fit to data 
in the signal region, as shown in Fig. [I] For this combined 
fit, the simulation and instrumental background shapes 
are fixed to those from previous fits, and only the abso- 
lute scale of the two types of background is allowed to 
float. The normalization of the background for simulated 
(MC) processes is found to be 1.06 ± 0.02 (statistical er- 
ror), in good agreement with the expected cross sections. 
The invariant mass distribution of the two leading jets 
after final background normalization is shown in Fig. [2l 

The standard DO neural network b tagging algorithm 
employs lifetime based information involving track im- 
pact parameters and secondary vertices [lOj. We opti- 
mize the choice of b tagging operating points for best 
signal significance and require one tight 6-tag (b-t&g effi- 
ciency ~50% for a mistag rate of ~0.4%) and one loose 
&-tag (6-tag efficiency ~70% for a mistag rate of ~4.5%). 
Table U shows the number of expected events from MC 
and instrumental backgrounds along with the number of 
events observed in data, before and after b tagging. Af- 
ter b tagging, 134 ± 18 events are expected and 140 are 



>1400r 




Invariant Mass (GeV) 

FIG. 2: Invariant mass distribution of the two leading jets 
before b tagging requirements. 



TABLE I: Number of events after selections. 



Sample 


No 6-tag 


Double b-tag 


ZH(M H = 115 GeV) 


2.46 ± 0.34 


0.88 ±0.12 


WH{M H = 115 GeV) 


1.75 ± 0.25 


0.61 ±0.08 


Wjj 


5180 ± 670 


7.6 ± 1.4 


Wbb 


397 ± 52 


35.4 ±7.1 


Wee 


1170 ± 150 


9.3 ± 1.9 


Z(-> r + T-)jj 


107 ± 14 


0.25 ±0.05 


Z(~* vu)jj 


2130 ± 280 


0.63 ±0.12 


Z(-> T + T-)bb 


6.39 ± 0.83 


0.63 ±0.13 


Z(-^ vv)bb 


229 ± 30 


24.9 ± 5.0 


Z(-^> T + T~)CC 


12.8 ± 1.7 


0.18 ±0.04 


Z{— > vv)cc 


467 ±61 


4.9 ± 1.0 


tt 


172 ± 34 


29.1 ±6.1 


Diboson 


228 ± 25 


3.84 ± 0.50 


Total MC Bkg 


10100 ± 800 


117 ± 17 


Instrumental Bkg 


2560 ± 330 


17.2 ±3.4 


Total Bkg 


12700 ± 900 


134 ± 18 


Observed Events 


12500 


140 



observed. 

Further signal-to-background discrimination is 
achieved by combining several kinematic variables 
using a NN. Independent MC samples are used for NN 
training, NN testing and limit setting. The instrumental 
background contribution is not taken into account during 
training, as its inclusion does not improve the expected 
sensitivity. The signal sample used for training is a 
combination of the ZH and WH contributions. Events 
are weighted such that the total contribution from each 
sample is that expected after b tagging. The NN input 
variables are the invariant mass of the two leading jets 
in the event, A1Z between the two jets, pr of the leading 
jet, pt of the next-to-leading jet, $r, Tfi-T and Ht- The 
input variables are selected for their ability to separate 
signal and background and to provide good modeling of 
data. The NN outputs for signal, background and data 
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t Data 

I Instrumental Bkg 

I Diboson 
I I W+jets 
I | Z+jets 




°0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

NN output 

FIG. 3: NN output distributions for M H = 115 GeV after 
b tagging. The MC expectation for the Higgs signal is scaled 
up by a factor of 50. 



are shown in Fig. [3J 

Systematic uncertainties affect the expected number of 
signal and background events ( "overall uncertainties" ) as 
well as the shape of the distribution in the NN output 
("differential uncertainties"). We estimate overall sys- 
tematic uncertainties associated with luminosity (6.1%), 
trigger efficiencies (5%), jet identification (5%), b tag- 
ging (7%), background MC cross section (6-18%) and in- 
strumental background (20%). All systematic uncertain- 
ties are common and correlated between signal and back- 
grounds, except for the uncertainties on the cross sections 
and the instrumental background. Differential uncertain- 
ties are estimated from the difference in the shape of the 
NN output by varying the jet energy scale (JES) by its 
uncertainties in a correlated way for all signal and back- 
ground MC samples at each mass point. The difference in 
the distribution of the NN output from the uncertainty in 
the shape of the MC di-6-jet mass spectrum is also taken 
into account at each Mh point. The JES uncertainty was 
estimated to be < 10% and that for the mass spectrum 
< 8%. Additionally, the impact on the NN output of the 
possible discrepancy in the low mass region in Fig. [5] was 
investigated and found to be negligible. 

We set a limit on the Higgs production cross section 
using a modified frequentist approach with a Poisson 
log-likelihood ratio (LLR) statistic 
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distribution is used to construct the LLR test statistic. 
The impact of systematic uncertainties is incorporated 
through "marginalization" of the Poisson probability dis- 
tributions for signal and background, assuming Gaussian 
distributions. We adjust each component of systematic 
uncertainty by introducing nuisance multipliers for each 
and maximizing the likelihood for the agreement between 
prediction and data with respect to the nuisance param- 
eters, constrained by the prior Gaussian uncertainties for 
each. All correlations in the systematics are maintained 



TABLE II: Expected (Exp.) and observed (Obs.) limits in pb 
and as a ratio to the SM Higgs cross section (in parentheses) , 
assuming H — * bb. 



Higgs Mass (GeV) 


105 


115 


125 


135 


ZH Exp. 
ZH Obs. 


1.6 (15) 
1.5 (14) 


1.5 (19) 
1.5 (20) 


1.4 (29) 
1.4 (30) 


1.2 (47) 

1.3 (51) 


WH Exp. 
WH Obs. 


4.8 (25) 
4.4 (23) 


4.3 (33) 
5.0 (39) 


3.8 (47) 
4.4 (55) 


3.6 (84) 
4.2 (99) 


VH Exp. 
VH Obs. 


2.8 (9.1) 
2.6 (8.7) 


2.5 (12) 
2.7 (13) 


2.3 (18) 
2.5 (20) 


2.0 (30) 
2.3 (34) 
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Expected Limit 

Observed Limit 

— SM Prediction 



110 120 



130 140 
Higgs Mass (GeV) 



F_IG. 4: 95 % C.L. upper limit on a(pp VH) x B(H -> 
bb) (and corresponding expected limit) for VH production vs. 
Higgs mass. 



between signal and background. The resulting limits are 
presented in Tabic HTl 

In summary, we have performed a search for the stan- 
dard model Higgs produced in association with either a 
Z or W boson (denoted as VH), in the final state topol- 
ogy requiring missing transverse momentum and two b- 
tagged jets in 0.93 fb -1 of data. In the absence of a 
significant excess in data above background expectation, 
we set limits on cr(pp — > VH) x B(H — > bb) at the 95% 
confidence level of 2.6 pb - 2.3 pb for Higgs boson masses 
in the range 105 - 135 GeV. The corresponding expected 
limits range from 2.8 pb - 2.0 pb. The expected and ob- 
served limits, along with the SM prediction, are shown 
in Fig. [4] as a function of Higgs mass. This is the most 
stringent limit to date in this channel at a hadron col- 
lider. 
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